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Changes in atmospheric circulation over the past five decades 
have enhanced the wind-driven inflow of warm ocean water 
onto the Antarctic continental shelf, where it melts ice shelves 
from below 1-3 . Atmospheric circulation changes have also 
caused rapid warming 4 over the West Antarctic Ice Sheet, 
and contributed to declining sea-ice cover in the adjacent 
Amundsen-Bellingshausen seas 5 . It is unknown whether these 
changes are part of a longer-term trend. Here, we use water- 
isotope (8 1s O) data from an array of ice-core records to place 
recent West Antarctic climate changes in the context of the 
past two millennia. We find that the 8 18 0 of West Antarctic 
precipitation has increased significantly in the past 50 years, 
in parallel with the trend in temperature, and was probably 
more elevated during the 1990s than at any other time during 
the past 200 years. However, 8 18 0 anomalies comparable to 
those of recent decades occur about 1% of the time over 
the past 2,000 years. General circulation model simulations 
suggest that recent trends in 5 18 0 and climate in West 
Antarctica cannot be distinguished from decadal variability 
that originates in the tropics. We conclude that the uncertain 
trajectory of tropical climate variability represents a significant 
source of uncertainty in projections of West Antarctic climate 
and ice-sheet change. 

The West Antarctic Ice Sheet (WAIS), which is grounded largely 
below sea level, is potentially unstable. Mass loss from the WAIS 
is contributing to present sea-level rise, owing to the widespread 
thinning of ice shelves and the acceleration of the large outlet 
glaciers that drain the ice sheet into the ocean 1 . Contemporaneous 
with the loss of mass from the WAIS, air temperatures over the 
WAIS have increased significantly in the past 50 years 4,6,7 . 

Climate and ice-sheet changes in West Antarctica are closely 
linked with one another by changes in regional atmospheric 
circulation 8 . Observations beneath the floating ice shelf of Pine 
Island Glacier, a major drainage system for the flow of the WAIS 
into the Amundsen Sea, show that the primary cause of ice-shelf 
thinning is the presence of warm Circumpolar Deep Water on 



Figure 1 1 Map of West Antarctica. Ice-core locations are shown by filled 
circles. Blue shading shows the main Siple Coast and Amundsen Sea ice 
streams. Ice shelves are shaded grey. The inset map shows the locations of 
well-dated, annually resolved ice cores for which there are & 18 0 data 
available to at least 1994. WD, WAIS Divide ice core (white-edged circle). 
PIG, Pine Island Glacier. The location of the Byrd weather station is shown 
by an open circle. 

the Antarctic continental shelf 3 . Circumpolar Deep Water inflow 
onto the continental shelf probably increased between the 1980s 
and 1990s because of increased wind stress at the shelf edge 2,8 . The 
patterns of sea-level pressure and geopotential height anomalies 
associated with increased westerly wind stress 8 also favour reduced 
sea-ice extent 9 and the advection ofwarm air onto the continent 4,10 . 

It is unknown whether the climate and glaciological changes 
that have occurred in West Antarctica in recent decades are part 
of a longer-term trend associated with anthropogenic climate 
forcing. This question cannot be evaluated with direct observations. 
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Figure 2 | West Antarctic temperature and & 18 0. a, Temperature 7 at Byrd 
Station (80° S, 120° W) (upper) and composite of annual mean & ls O 
anomalies from ice cores in West Antarctica (lower). The two time series 
are correlated at r= 0.48, p = 0.01. Grey shading shows the running 
decadal mean of temperature and the standard error of the running decadal 
mean of & 18 0. b, Probability from a one-tailed f-test that decadal mean 
West Antarctic S 18 0 centred on any given year is as elevated as the decade 
of the 1990s (1991-2000). The dashed line shows the number of records 
contributing to each decadal mean. 


West Antarctic temperature and pressure observations begin only 
in 1957, and reliable satellite observations of Antarctic sea ice 
date to 1979. Comprehensive observations of glacier dynamics 
in the most rapidly changing areas were initiated in the 1990s. 
Borehole temperature data from the WAIS, although confirming 
the recent rapid rise in temperature, do not resolve decadal-scale 
variability in the past 6 . 

Here, we use records of the oxygen isotopic composition 
of precipitation (8 ls O, see Methods) from ice cores to assess 
the significance of recent West Antarctic climate trends, and to 
place the observed glaciological changes in a longer-term context. 
Unlike temperature or other conventional climate variables, S ls O is 
relatively well sampled in West Antarctica and thus provides the best 
available data for this purpose. We use a network of ice cores (Fig. 1 ) 
that includes sixteen new and updated multi-decade to century- 
length records from across the WAIS, and a new 2000-year-long 
record from the central ice divide (WAIS Divide). 

The use of S ls O in precipitation as a proxy for temperature 
is well known, and supported by our data. However, we do not 
use 8 ls O as a proxy for temperature. Rather, we observe that 
S ls O in West Antarctica covaries with atmospheric circulation in 
a manner similar to temperature. Positive anomalies in both near- 
surface temperature and S ls O of precipitation in West Antarctica 
arise when anomalous meridional flow in the troposphere advects 
warmer air from lower latitudes, reducing the net distillation of 
water vapour that otherwise progressively lowers 6 ls O values in 
precipitation 11 . The 8 ls O of precipitation in West Antarctica is also 
increased by reduced sea-ice cover, owing both to an increase in the 
fraction of locally derived moisture, and an increase in boundary- 
layer turbulence that produces more effective penetration of local 



Figure 3 | Decade-average 6 ls O from the WAIS Divide ice core for the 
past 2,000 years. Grey shading shows 2 s.d. about the decadal mean, 
based on the upper 100 years of the multi-core & 18 0 composite, providing 
an estimate of the 95% confidence range. The dashed line shows the 97.5 
percentile value relative to the average linear trend. 

maritime air masses onto the continent 12 . West Antarctic 8 ls O thus 
reflects the same circulation anomalies that have contributed to 
the trends in temperature, sea ice and the glaciologically significant 
ocean circulation changes of the past few decades. 

Figure 2 shows 6 ls O averaged over the West Antarctic ice-core 
records for the past 200 years. A prominent feature is the elevated 
5 18 0 of the 1990s and significant upward trend since the 1950s, 
consistent with the surface temperature trend 4,6,7 . In nearly all 
of the individual ice-core records, the highest 8 ls O values occur 
after 1985, and the mean decadal maximum occurs during the 
1990s (1991-2000). The limited available data indicate lower values 
during the 2000s. The mean upward trend in 6 ls O over the past 50 
years therefore primarily reflects the anomalous values ofthel990s. 

Before the 1990s, there are two other prominent decadal-scale 
anomalies in West Antarctic 8 ls O in the past two centuries: the 
1830s and the 1940s (1936-1945). We use the population statistics 
from the ice-core records to assess the level of confidence at which 
the S ls O values of these decades, averaged over West Antarctica, are 
distinguishable from the 1990s maximum (Methods). The results 
(Fig. 2) show that decadal-average 8 18 0 since 1990 is probably more 
elevated than that of the 1830s and 1940s, but only at low confidence 
( p ~ 0.2 and p ~ 0.3, respectively). It is extremely likely ( p < 0.05) 
that the mean 8 ls O of the decade 1991-2000 was more elevated 
than during any decade in the previous 40 years, consistent with 
the evidence for widespread warming and large-scale changes in 
atmospheric circulation 4,10 . 

Analysis of the long record at WAIS Divide shows that 5 18 0 in 
West Antarctica is anomalous not only with respect to the past two 
centuries, but also with respect to the past two millennia. WAIS 
Divide is well situated to be representative of West Antarctica as 
a whole: the WAIS Divide 8 ls O record is equally well correlated 
with 6 ls O records from cores east and west of the ice divide 11 . 
Decadal-average S ls O values comparable to the 1990s in the WAIS 
Divide record are reached on only four occasions in the past 1,000 
years (Fig. 3). Before 1,000 years ago, modern decadal-average 6 ls O 
values are reached more frequently, but these are superimposed on a 
declining trend attributable to the influence of Milankovitch orbital 
forcing and ice flow 13 . Assuming that the decadal variability is 
independent of orbital forcing, we calculate 6 18 0 anomalies relative 
to the long-term trend (dashed line in Fig. 3). Anomalies in 6 ls O 
similar to those of the 1990s occur just twice in the past 2,000 
years; assuming sampling error estimated from the multiple shorter 
records, comparably elevated S ls O values were reached about 1% 
of the time. Analysis of trends yields a similar result: trends in 
5 18 0 as large as the most recent 50-year trend in the WAIS Divide 
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Figure 4 | Modelled versus observed West Antarctic 6 ls O and tropical 
SSTs. a, Difference in mean simulated decadal-mean & 18 0 in Antarctic 
precipitation between 1991 and 2000 and the three preceding decades, 
b, Model range (2 s.d., grey shading) of simulated annual mean 6 18 0 in 
precipitation averaged over West Antarctica compared with observed & ls O 
anomalies (black line). In both a and b, the model simulations are from 
10-member ensembles of ECHAM4.6 simulations forced by global tropical 
SST with a slab ocean in the extratropics, c, SST anomalies 26 averaged over 
the central tropical Pacific Niho4 region (thin solid line) and over the entire 
tropics (thick line), 1880-2009. 

record is reached ~2% of the time in the past 2,000 years. Hence, 
the anomalous S 18 0 values of recent decades in West Antarctica 
do not seem to be unprecedented, but are near the upper limit 
of the range of natural variability. A similar conclusion has been 
reached for the northern Antarctic Peninsula, based on data from 
an ice core there 14 . 

Recent climate trends in West Antarctica have been shown to 
result from anomalous convection over the tropical Pacific, propa- 
gated to the high latitudes by atmospheric Rossby waves 9,10,15 . Sig- 
nificant changes have also occurred in the large-scale circumpolar 


westerlies, expressed as an increase in the Southern Annular Mode 
(SAM) index, in austral summer 16 . Changes in the SAM may be 
reflected in dust and aerosol increases recorded in West Antarctic 
ice cores 17 . However, the largest trends in both temperature and 
in ice-core 8 ls O in West Antarctica have occurred in the non- 
summer seasons, during which the influence of tropical forcing 
has dominated 15 . Comparisons of West Antarctic 8 ls O records 
with tropical sea surface temperature (SST) and rainfall records 
suggest an important role for the tropical Pacific in influencing West 
Antarctic decadal 8 ls O variability throughout the past 200 years 18,19 . 
Indeed, the statistics of climate variability in West Antarctic inferred 
from 8 ls O are similar to those for tropical Pacific climate. Although 
the 1 990s were the warmest decade of the past century in the tropics, 
they are distinct from other decades only at moderate confidence 20 
and proxy SST data indicate that conditions similar to those of the 
1990s in the tropical Pacific occur with a frequency of about twice 
per century over the past 350 years 21 . 

The statistics of the ice-core 8 ls O records suggest that decadal 
variability in the tropical Pacific is sufficient to account both for the 
strong recent trends, and for the magnitude of decadal variability 
in West Antarctic climate over the past century. This is supported 
by simulations with a general circulation model (GCM). We used 
the ECHAM4.6 GCM (ref. 22), which has been shown previously 
to reproduce the magnitude and spatial pattern of observed 
atmospheric circulation and temperature change in West Antarctica 
in response to tropical SST boundary conditions 10 . When run 
with its water-isotope module 23 , ECHAM4.6 reproduces the 
characteristic spatial covariance relationships between 8 18 O in West 
Antarctic precipitation and geopotential height anomalies and SST 
seen in observations (Supplementary Information). When forced 
either by global observed SSTs or tropics-only SSTs, ECHAM4.6 
captures the magnitude and spatial expression of the prominent 
1990s decadal 8 ls O anomaly (Fig. 4). As in the observations, the 
simulated 8 ls O of the 1990s is significantly elevated compared 
with the past several decades, but is only marginally distinct from 
other decadal maxima over the past century. Furthermore, whereas 
there is a significant positive trend in the tropical SST boundary 
conditions in the past century, there is no significant trend in 
the resulting S ls O simulated over West Antarctica, consistent with 
the observations. This occurs because it is the pattern of warming 
in the tropics, rather than the mean temperature, that produces 
the Rossby-wave response 10,24 . The correlation between the forcing 
(tropical SST) and simulated 8 ls O time series is only marginally 
significant {r = 0.22, p = 0.07), whereas the correlation between the 
central tropical Pacific SSTs and simulated 8 ls O is highly significant 
(r = 0.58, p < 0.001). 

Much attention has been given to the role of anthropogenic 
radiative forcing on Antarctic climate through the influence of both 
increasing greenhouse gases and stratospheric ozone decline on 
the SAM (ref. 16). Our results suggest that in West Antarctica the 
influence of natural decadal variability originating in the tropical 
Pacific has masked any anthropogenically forced response. A caveat 
is that there is a small but robust long-term trend in the pattern 
of tropical SSTs that may be interpreted as a response to radiative 
forcing 25,26 . However, most analyses suggest that observed changes 
in the key central Pacific and El Nino/Southern Oscillation regions 
that influence West Antarctic climate cannot be distinguished from 
stochastically forced variability 20,26 . 

Our results have implications for the response of the West 
Antarctic Ice Sheet to anthropogenic forcing. How tropical Pacific 
climate will change in the future is uncertain, with some GCMs 
projecting a more El Nino-like state and others suggesting a 
tendency towards greater La Nina-like conditions 27 . These different 
possibilities have potentially quite different impacts on the regional 
circulation affecting West Antarctica, suggesting that tropical 
climate represents a significant source of uncertainty in West 
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Antarctic climate, and in the contribution of the WAIS to sea level, 
over the next century. Moreover, some GCMs overestimate the 
magnitude of variability associated with the SAM (ref. 28) — the 
dominant response to radiative forcing in most models — and are 
therefore likely to underestimate the importance of variability 
associated with changes in the tropical Pacific. Indeed, large- 
ensemble modelling studies show that internal variability associated 
with the tropics is a significant source of uncertainty in climate- 
change projections at high latitudes, especially for the atmospheric 
circulation 29 . Projections of the contribution of the WAIS to future 
sea-level rise that are based on present rates of ice-sheet mass loss 30 
should be treated with caution. 

Methods 

Ice-core data presented in this manuscript include both previously published 
and new analyses of ice and firn cores, primarily from the United States’ ITASE 
(International Trans- Antarctic Expedition). The development of depth-age scales 
for these cores was conducted using multi-parameter high-resolution chemistry. 
Water stable-isotope concentrations were measured by mass spectrometry and 
laser spectroscopy at the University of Washington, and are reported as 6 ls O 
deviations in per mil (%o) from Vienna Standard Mean Ocean Water (VSMOW), 
normalized to VSMOW and SLAP (Standard Light Antarctic Precipitation), 
where 6 ls O = R/Rvsmow — 1, R is the abundance ratio of 18 0/ 16 0 in water, 
6 18 Ovsmow = 0.0%o and 6 18 Oslap = — 55.5%o. The ice-core 6 ls O composite records 
were produced by normalizing the mean over the data overlap period common to 
all cores. The spatial distribution of decadal-mean 6 ls O data from West Antarctic 
cores is indistinguishable from a normal distribution (p < 0.05) based on a Lilliefors 
goodness-of-fit test. We use a one-tailed Student’s f-test of the null hypothesis that 
a given decade mean comes from the same population as the 1990s (1991-2000) 
mean. Use of the non-parameteric Wilcoxon rank-sum test yields indistinguishable 
results. The significance level p of reported correlations accounts for the degrees of 
freedom reduced by the quotient ( 1 — r\ r 2 )/( 1 + r\ r 2 ), where r,- is the lag- 1 autocor- 
relation of time series i. The general circulation modelling results are from exper- 
iments with the ECHAM4.6 atmospheric GCM, at a horizontal resolution of T42 
(~2.8° latitude x 2.8° longitude) with 19 vertical levels, using observed SST (ref. 26) 
as a boundary condition. In experiments using only tropical SSTs as the boundary 
condition, the atmospheric GCM is coupled to a slab ocean in the extratropics. 

Data, http://nsidc.org/data/nsidc-0536.html. 
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